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A non linear model associated with a Landau-Ginzburg-like behavior in mean field 
approximation forecasts phase transition waves and solitary kinks near the critical 
point. The behavior of isothermal waves is different of the one of isentropic waves 
as well in conservative cases as in dissipative cases. 



1. Introduction 

The study of plane waves propagating in fluid second gradient theory is 
a main subject in continuum mechanics From our paper (®) where 

strong variation of density for the matter constituent occurs in isothermal 
conservative motions, we propose an extension to other thcrmodynamical 
conditions and for dissipative motions. Such cases appear in phase transi- 
tions through interfacial layers Near the critical point the thickness 
of the interface separating the two phases of a fluid gets non molecular di- 
mension. The theory of second gradient corresponding to an extended van 
der Waals model with thcrmodynamical potentials modeling the mean field 
approximation allows to study the interface in equilibrium Wi^jj a 

convenient rescalling, the governing equations of material waves can be ex- 
pressed for conservative and viscous motions with a Landau-Ginzburg-like 
behaviour they are compatible with the second law of thermodynam- 
ics and we can rewrite the one-dimensional equations of non homogeneous 
fluids in terms of density i^^^^. To obtain wave motions, we study the 
thcrmodynamical potentials of a fluid near the critical point and extend 
the expression given in ^'^^ to the local internal energy. Thcrmodynamical 
potentials are introduced in equation of motions in two cases: when the 
temperature is fixed (isothermal motions) and when the specific entropy 
along trajectories is constant (isentropic motions). We find the possibility 
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of non-linear waves in form of solitary waves or phase transition waves for 
conservative motions, but no kinks as traveling waves can appear for non 
conservative motions. Finally we discuss the fluid velocity near the critical 
point depending on the boundary conditions. 

2. Thermodynamic potentials near the critical point 

In mean field approximation, the thermodynamical potentials of fiuids can 
be expressed in an analytical form near the critical point The chemical 
potential of a homogeneous fluid can be written in the following form ^'^^ 

li = ll' + tll^ (T - Te) + (p - p,) (T - Te) + i {P-Pcf (1) 

where p — Pc and T — T^ denote the difference between the density p and the 
Kelvin temperature T from their values Pc, at the critical point. This 
simple expansion depends on coefficients corresponding to the values of the 
partial derivatives of p, at the critical point; that is to say, 

d^+i p 

The chemical potential p has been expanded using a cubic polynomial in 
terms of density and temperature; they are not the natural variables but 
the most convenient for the calculations. ^ 

Denoting the coefficients hy A= pl-^^ > 0, B = - p%Q > (the sign condi- 

6 

tions are at the critical point) and p,o = p," + /Uqi {T — Tc) depending only 
on T, Eq. (1) reads 

M = Mo - A{n -T){p- p^)+B{p- p^f (2) 

Expression (2) allows us to solve the problems of waves in an universal way 
where T—Tc plays the role of an order parameter. Such a chemical potential 
expansion agrees with the van der Waals equation of state where only 
the values A and B are representative of a particular ffuid. 
Consequently, the free energy per unit volume ip can be written in the form 

= MoP - ^{Tc -T){p- p,f + |(p - - /(T) (3) 

where / is an additive function of T only. The pressure p verifies p = pp—tp 
and we get near pc the expansion 

p = Pc + DiT-Tc)-iv-Vc) (Ai (T - Te) + Biiv- v^)) (4) 
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1 dj) 
where v = — denotes the specific volume, D = {vc,Tc), f{Tc) = pc 

and Ai = A pi, Bi = B pi. 

Due to the general properties of two-phase regions near the critical point 
the two-density version of specific internal energy of the fluid is expressible 
in the form 

a{v,s) = {x^ -yf + + Px + -/y + 5 (5) 

with X = ar + ba, y = cr + da where r = v — Vc, a = s — Sc , 
7, 6, a, b, c, d are constants and s denotes the specific entropy of the fluid. 
Takmg mto account that p = — —(v,s), J = -^(v,s) and denotmg 

OV OS 

ai{v, s) = a{v, s) + Pc{v — Vc) — Tc{s — Sc) — r^, we obtain the partial 

differential equation 

dai ,^ , , 9ai 
^ + iD-A,r)—=0 

which yields the general solution for the specific internal energy 

Ai 

a = g{a - Dr + — r'^) + Bi— - p^r + T^a 

Due to the fact that a must verify the general form (5), g{u) is necessarily a 

Ai 

polynomial with respect to u = a—Dr+-^ in the form g{u) = Eu+Gu^. 

The term Gu^ yields a term of four degree in r and r must be a function of 
X only. No term in the form xy appears in (5); consequently a — Dr does 
not contain the variable y and the general form of a is 

a = E{k2 y+^klx'') + G{k2 y+^kl x'^f + ^ktx^ - p^r + T^(j 
where E and G are two constants. 

By straightforward calculations, an identification with the expression (5) 
yields the general expression for the specific internal energy 

a{v, s)=x^- 2x'^y + 2y'^ - p^{v - v^) + T^{s - s^) (6) 

with X = \^ {v- Vc) and y = -^J ^ [D{v - Vc) - {s - Sc)] 
V 2 A.I V 2 

As a consequence the pressure gets the form 

/ B-\ , 0^ D / . O 



P=Pc + 4{/ -^{y-x')x + ^ V2Bi K - %) (7) 
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and in the case where s is constant, the specific enthalpy reduces to the 
form 

where ho is a constant depending on s. 

Near the critical point, a fluid behaves like a gas but with a high density 
As for interfaces separating two bulks, the fluid is not homogeneous. The 
view that a non-homogeneous fluid near its critical point may be treated 
as matter in bulk with a local energy density that is that of a hypotheti- 
cally uniform fluid of composition equal to the local composition with an 
additional term arising from the non-uniformity, and that the latter may be 
approximated by a gradient expansion typically truncated in second order 
is most likely to be successful and perhaps even qualitatively accurate i^'i^. 
The first study has been done on the theory of the near-critical interface 
within the framework of the van der Waals theory of capillarity ^; the sim- 
plest model able to take into account the density and its gradient uses a 
unique supplementary quantity represented by the constant C of internal 
capillarity. In S.I. units, the value of C for water at 20° Celsius is of the 
order of 10~^^ (see ^^). In the mean field approximation the expression of 
the internal energy anh of a non homogeneous fluid near the critical point 
is in the form^O'i^." 

C 

anh = a+ — {Wpf 

The supplementary term due to the non-homogeneity of the medium in the 
expression of the internal energy is effective only in interfaces and for fluids 
near the critical point. 

3. Motions of a fluid near the critical point 

Equations of motions are classically given in the literature with the additive 

second gradient term, Fluid motions are associated in the literature with 
the capillary fluid equation of motion 

pr + Vp + pV(0-CVV)-divcr^, =0 (9) 

where p is the previous thermodynamic pressure associated with the 
medium considered as homogeneous, (i.e. the pressure entering in the ther- 
modynamic potentials of the homogeneous fluid near the critical point), Q. 
is the body force potential, (T„ = Atr(A) Jd-I- 2/iA is the viscous stress 
tensor in the classical form where A is the velocity deformation tensor. 
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r the acceleration vector, V the gradient operator and the Laplacian 
operator. The mass balance yields 

|^+divpV = (10) 
where V denotes the velocity of the fluid. 

Let us notice that we can obtain the equation of energy such that the sys- 
tem is compatible with the second law of thermodynamics i^-^o. Capillary 
fluids belong to the class of dispersive systems because the internal energy 
depends not only on density but also on its derivatives with respect to space 
variables; we have previously seen that the constant solutions are neverthe- 
less stable 

Let us study the one- dimensional problem when the velocity V and the 

density p arc only functions of the variable C = z — ci, where z is the space 
variable, t the time and c the wave celerity with respect to a Galilean frame: 

V = V{z - ct), p = p{z-ct) 

The mass balance equation yields 

, d{pV)_ 

and by integrating, we obtain 

p{V-c) = q (11) 
where q is constant in the motion. In the cases of waves we obtain 

We consider that body forces are negligible as in space with weightlessness. 
In the uni-dimensional case, 

In the following, we assume that v = (A + 2 /u)/p is constant in the fluid as 
assumed in interfaces^''. Then Eq. (9) yields 

We consider two cases a) isothermal motions corresponding to fluid motions 
near equilibrium conditions and b) isentropic motions corresponding to fast 

1 . . T 1 . 1 dp dH , „ . ^ 

velocities. In the two cases, we can write -— = — — where H is the 

p dC d( 



div ay = {X + 2p)q^^^ 
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chemical potential of the fluid in case a) and the specific eiitlialpy in c;asc h). 
In all the cases, near the critical point p w Pc, consequently vq/ p^ ~ fq/p^ 
and the equation of motion becomes 

c'^.-^'i=H-H. + \'-, (12) 



where Ho is constant. 



4. Example of waves of a fluid near the critical point 

Two main cases arc generally considered for traveling waves : isothermal 
processes and isentropic processes. Motions can be conservative as a math- 
ematical limit of the dissipative case. Nevertheless, in thin interfaces the 
viscosity may be neglected and conservative cases may be considered as 
realistic physically . The cases of viscous fluid involve inequality due to 
Liapounov functions. 



4.1. Existence of solitary waves 

Let us notice that for viscous capillary fluid, solitary waves as kinks cannot 
appear. Multiplying Eq. (12) by ^ we obtain by integration : 

where Co is a constant and K'{p) = H{p) — Hg. 

Eq. (13) allows to obtain a first integral only when v = Q . Due to the fact 
that 

for z/ > and C > Co , j^^ y > 

it is not possible to obtain p(— oo) = p(+oo) and no solitary wave can 
appear in dissipative motions. 



4.2. Isothermal waves 

In the isothermal case the wave motion is obtained by using the chemical 
potential p, given by Eq. (2) 

1 (P n V q dp A ,^ „^ , ^ B , q^ , , ^ 

2 d^-2^ i - -2C ^^-^^ (^-^^^+2C (^-^^)'+4^ (1^) 
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where ko is a constant. 

Near the critical point \{p— pc)/ Pc\ ^ 1; consequently, we get the following 
expansion to the second order in {p — Pc)/Pc 



11/ 2{p-p,) , 3(p-pe) 



p^ Pi V Pc Pi 



and to the fourth order Eq. (14) yields 

1 d^p vq dp B / q^ 



2dC 2piCdC 2Cy 2Bpi 

where ki is constant. We define the following change of variables: 

^ ' (l+gF), C = i^eand Q= (15) 



2BpiJ' ^ ^ ^ (2C7)3/2Bi/2p2 

2C ^(^c-T) + ^ 

with ^ ill and ^ Pc 



Pt V 2Bp4 

Then, the equation of isothermal waves writes 
1 d^V dV 

Conservative case, y = Q: we obtain the first integral 

^)' = (l-FY-air-6i (17) 

where Oi and bi are two constants. The intersections of the straight line 
aiY + bi and of the quartic (1 — y^^a yjgj^j ^j^g density range 

Interfacial propagation - We recall the main results of (*): in the bulks 
phases, the densities are constant and thus the first and second derivatives 
of Y are zero. The straight line is tangent to the quartic at the associated 
points. If the bulk phases are different on both sides of the interface, the 
straight line has to be bitangcntial to the quartic which implies oi and bi 
are null (case of two phases). A liquid- vapor interface wave is similar to the 
one obtained in the equilibrium case {q = 0) but mass flows through the 
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interface. Vaporization or condensation phenomena depend on the sign of 
q. This case corresponds to a shock wave in the sense of Slemrod ^ and 

(see left side of fig. 1). 

Solitary waves - It is also possible to obtain a solitary wave moving in 
one of the bulk phases (liquid or vapor). The straight line is tangent to 
the quartic at the point associated with the bulk Y = Yo and intersect the 
quartic at the point associated with the middle of the wave Y = Y^. When 
Y belongs to the interval between Yg and Y^, the right hand side of Eq. 
(17) must be positive and this authorizes only two possibilities: 
a) — 1 < Yo < — 1/-\/3 and Y^ <Yl^ < 1, the density increases in the wave 
(see right side of fig. 1), 

h) l/\/3 < < 1 and —1 < Y^ < Yg, the density decreases in the wave. 




5 



Figure 1. On the left side of the figure a phase transition wave is represented (densities 
pi and pg are associated with the liquid and vapor bulks); on the right side a traveling 
wave is a kink associated with case a) (density po corresponds to the bulk of the fluid). 



Consequences - The change of variables (15) assumes that Tc — T+ — — 

is positive; Then, for T > T^, the flow through the interface must be more 
important than the limit qm = sj Ap'l{T — Tc) and the wave velocity £ must 
be greater than the limiting value ^rn = \J Apc{T — Tc) which is the celerity 
of isothermal acoustic waves. 

Dissipative case, > : In this case no solitary wave appears (see 
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section 4.1). Let us define Z{Y} =^ ——; then Eq. (16) yields 

^ZZ' -vQZ = Y^ -Y + kx 
We look for the phase transition waves with a solution in the form 

Z{Y) = aiY^ +(.hY + 

where the polynomial aiY'^ + (iiY + 71 has two real roots. Then, the 
differential equation 

dY 

— = aiF^ + + 71 

has a solution in the same form than (18); (see left side of fig. 1). 
Straightforward calculations prove that the polynomial has two real roots 
when 

(^A (Te - T) + < 24 C^Bpl 

which corresponds to a velocity q/ pc of the wave small enough. This con- 
dition does not appear in the conservative case where i/ = 0. 



4.3. Isentropic waves 

For an isentropic conservative motion, s = Cte and the equation of motion 
is in a form deduced from Eq. (12), 

where h is now the specific enthalpy. We consider the new change of vari- 
ables 

p = Pc{l + Y), ^ = LC and q = bQ 

with 

DBi ' pcAi 
and the equation of wave motions yields 

i ( d^j =^(l + 3Q') + ^(r-2Q2) + aoF + 6„ (19) 
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2pcD 

with T = — - — and ao and ho are two constants. We denote by = 
^1 

-"2 — Wc a limit celerity of the waves with respect to the fluid. Due to 

the form of the second member of Eq. (19), we notice immediately that 
it is not possible to obtain phase transition waves. Consequently traveling 
waves cannot appear in dissipative motions when the entropy is constant 
in the fluid. The two cases are represented on fig 2. 



f(p) f(P) 




Figure 2. On the left side of the figure, the case where the Une fiY) = cioY + bo is 
above the cubic representing the left side of Eq. (19) is presented; on the right side the 
opposite situation: the line f{Y) = aoY + bo is below the cubic. 



Case 1: < r/2 or i < im- In this case solitary waves cannot appear 
(see left side of fig 2). 

Case 2: > t/2 or £ > In this case solitary waves are possible 
depending on initial conditions (sec right side of fig 2). 
In the van der Waals model of pressure we obtain by straightforward cal- 
culations 

P, , QPcPc J, 3pcPc 



and consequently, £m = \ ■ 

\ 3 Pc 

4.4. Fluid velocity near the critical point 

For a given density, the fluid velocity is deduced at time t = from Eq. 
(11). Consequently, 
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where V{z) is the fluid velocity at z and c is an arbitrary constant. At any 
time t, 

V{z-ct) = c+ ^ 

p[z — ctj 

where initial conditions yield the arbitrary velocity c. 

Interface propagation - In a phase transition wave, the fluid changes 
from liquid to vapor as its volume increases. Consequently such a phe- 
nomenon cannot occur in a closed tube but only in a tube with only one 
closed end. For example, the tube is closed at the other end with a piston 
whose displacement is imposed. The fluid velocity at the fixed end is zero 
(for example liquid bulk) and c = —q/pi- If we impose a value U for the 
piston velocity in the vapor bulk, we deduce the value of q 




where and pi are the values of the density in the vapor and liquid bulks. 
From the values of p^ and pi deduced from expression (18), we determine 
the flow q and the velocities f ~ q/Pc and c. 

Solitary wave - The volume of the interface that moves in only one 
bulk phase is constant. Such a wave may be moving in a closed tube 
such as Natterer tube . At the ends (assumed far from the region of 
the wave) the velocity in the bulk phase is zero and c ~ —q/pc- When 
the temperature is close to Tc, the model of fluid endowed with internal 
capillarity allows to obtain traveling waves in a tube which depends on two 
arbitrary parameters. 

5. Conclusion 

The mean field approximation for fluid near the critical point is able to 
predict solitons and transition of phase waves. Two kinds of waves are 
investigated: a) liquid- vapor waves in the case of isothermal medium whose 
the celerity depends on the distance between the temperature and its critical 
value. Above the critical temperature, the fluid behaves like an elastic 
medium and the waves are supersonic with respect to the isothermal sound 
velocity ^. Below the critical temperature, the fluid behaves like a nonrigid 
medium and the wave velocities can take any value in conservative motion 
but are bounded in dissipative motion proportionally to the viscosity of 
the fluid, h) Traveling kinks appear in the conservative case as well as for 
isothermal than for isentropic motions. 
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Thermodynamics functions and parameters vanish or diverge at the critical 
point proportionally to some power of the distance from that point currently 
measured as T — Tg. The critical exponents are central to the discussion of 
critical phenomena and can be generalized by nonclassical value of critical 
point exponents for the potentials given in section 2. The problem can 
be extended for multi-component fluid mixtures but critical points are not 
unique; for a mixture of fluids there is a curve of critical points. 
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